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ABSTRACT 

Recently It has been argued based on theoretical calculations and 
experimental data that there Is a universal form for the equation of state of 
solids. This observation was restricted to the range of temperatures and 
pressures such that there are no phase transitions. The use of this universal 
relation to estimate pressure-volume relations (l.e., Isotherms) required 
three Input parameters at each fixed temperature. In this paper we show that 
for many solids the Input data needed to predict high temperature 
thermodynamical properties can be dramatically reduced. In particular, only 
four numbers are needed: (1) the zero pressure (P = 0) Isothermal bulk 
modulus; (2) Its P = 0 pressure derivative; (3) the P = 0 volume; and (4) the 
P = 0 thermal expansion; all evaluated at a single (reference) temperature. 
Explicit predictions are made for the high temperature Isotherms, the thermal 
expansion as a function of temperature, and the temperature variation of the 
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Isothermal bulk modulus and Its pressure derivative. These predictions are 
tested using experimental data for three representative solids: gold, sodium 

chloride and xenon. Good agreement between theory and experiment Is found. 

I. INTRODUCTION 

Knowledge of the equation of state (the pressure, volume, temperature 
(P-V-T) relation), based on either calculation or measurement, Is of primary 
Importance In both basic and applied sciences. It provides Insight Into the 
nature of solid state theories, and determines the values of fundamental 
thermodynamic parameters. 

Many theoretical and semi-empirical calculations have been carried out to 

describe the Isothermal equation of state (EOS). Analytic semi-empirical 

expressions for the EOS for each class of solid abound in a variety of forms, 

(Ref. 1), with care often taken to present different forms for different 

classes of solids (Refs. 1 to 3). We have found, however, that there is a 

universal EOS for all classes of solids In compression and In the absence of 

phase transitions (Ref. 4). Further, we found that this universal EOS can be 

predicted from a knowledge of three equilibrium quantities at each fixed 

temperature T: the equilibrium volume V , Isothermal bulk modulus B , 

o o 

and Isothermal (dB/aP) Q . Here, we use the notation (3B/3P) o for the zero 
pressure value of the Isothermal pressure derivative of the Isothermal bulk 
modulus: (dB/aP)^ p 

In the present study we will show how to predict high temperature 

properties of the EOS based on an even smaller amount of Input data: the 

thermal expansion coefficient a , V , B and (3B/3P) at zero 

ooo o 

pressure and at a single (reference) temperature. This approach, Involving 
Input data at a single temperature, applies to solids whose thermal pressure 
Is essentially Independent of the volume and linear with temperature above the 
Debye temperature. This property of near volume Independence and linearity 
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with temperature has been found for a number of solids (see, e.g., Refs. 5 
to 8). Based on this assumption the new approach not only predicts the 
pressure-volume relation at high temperatures, but also the temperature 
dependence of the thermal expansion, B q , and (3B/3P) q . We will see that 
rather good accuracy can be obtained with these predictions In comparison with 
experimental data for some typical solids (a metal, an alkali-halide and a 
rare-gas solid). Even nonlinear effects in the thermal expansion, which have 
been difficult to obtain (Ref. 9), are accurately predicted In this approach. 

In the following we will first briefly review the observation of a 
universal EOS for all classes of solids in compression (Ref. 4). Next we will 
formulate a method to predict high temperature properties of the universal EOS 
based only on zero pressure data at a reference temperature. We then present 
a simple analytical approximation to the universal EOS which is then used as 
the reference isotherm. The last section will show a set of predictions for 
the thermodynamic properties of some typical solids (a metal, an alkali-halide 
crystal, and a rare-gas solid). These predictions will Include high 
temperature pressure-volume curves, thermal expansion curves, and the 

temperature dependence of the Isothermal bulk modulus and Its pressure 

* . / 

derivative. 

II. A UNIVERSAL EQUATION OF STATE 
A universal energy relation between the total energy E and a 
characteristic distance has been discovered for bimetallic adhesion (Ref. 10), 
chemisorption on metals (Ref. 11) and metallic cohesion (Ref. 12) (and even 
then hypothesized for the EOS of nuclear matter (Ref. 13)). This relation can 
be written as 

E ( a ) = AE E*( a* ) ' (2.1) 
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where 

( r - r ) 

a* = v WS WSe' (2.2) 

ft 

Here AE Is' the equilibrium binding energy for metallic cohesion, r ws Is 
the Wlgner-Seltz radius with P = 0 value of r WSg and ft Is a scaling 
length. 

Consider now the classes of solids other than metals. There Is some 
evidence (Refs. 12 and 14) that covalent bonds satisfy Eq. (2.1). However, 
rare gas and Ionic bonds do not. For energy relations In either compression 
or expansion, Ionic and rare-gas Interactions are of Inherently different form 
than the wave-function-overlap Interactions (Ref. 15) characteristic of 
metallic and covalent bonds. It Is these differences that have led to the 
proposal of different forms for the EOS of different classes of solids 
(Refs. 1 and 2). 

However we have found that In compression the form of the pressure-volume 
relation Is dominated by overlap Interactions for all classes of solids. That 
Is, at T = 0 K, the pressure P(V) has the form 

. P( V) = - dE/dV (2.3) 

for all classes of solids, where E Is the universal form for the energy 

3 

(Eq. (2.1)) and V = 4wr ws /3 Is the volume. Thus many different forms 
are not needed to describe the EOS. Rather, a single universal expression 
suffices. 

In this paper we are primarily Interested In temperature effects. For 
T > 0 Eq. (2.3) becomes 

P(T, V) = - (dF(T,V)/aV) T = - dE( V)/dV + P*(T,V) (2.4) 

where F ( T , V ) Is the Helmholtz free energy, which Includes a temperature 
dependent Internal energy and entropy, and P*(T,V) Is the thermal pressure. 

Now by definition, P*(T,V) -* 0 as T •* 0. Thus for low T, - dE/dV 
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dominates and one might expect the form of P(T,V) to be similar to that of 
- dE(V)/dV with the temperature dependence being contained In the scaling 
parameters. For higher temperatures , P*(T,V) Is not negligible, but In 

general It Is known (Refs. 5 to 8, 16, and 17) to be a weak function of V. 

In this case, again one might expect that the form of - dE(V)/dV and P(T,V) 
would be similar with the scaling parameters containing the temperature 
Information. 

To make this more concrete, we first define the function H as 



This definition Is prompted by the form of - dE(V)/dV obtained from the 
universal energy relation, Eq. (2.1), for metals (Ref. 18). Our results for 
metals (Ref. 18) suggest that plots of In H versus 1 - (V/V Q ) 1//3 
should be approximately linear, of slope = 3/2 [(aB/3P) Q - 1] and Intercept 
= B q , where the subscript o refers to equilibrium values. We found 
(Ref. 4) that this Is In fact an accurate description of EOS data for many 
materials from all classes of solids over a broad range of temperatures. This 
Is an experimental confirmation of the expectations expressed above based on 
the weak volume dependence of P*(T,V). 

As further confirmation, we added a volume-independent P*(T) to 
pressure-volume data (Ref. 6) for cesium of strength sufficient to Increase or 
decrease V Q ( T ) by up to 4 percent. This thermal volume change and 
corresponding change In thermal pressure Is larger than one could actually 
encounter In solid cesium. Nevertheless, we found that the linearity of plots 
of In H(V,T) versus 1 - (V/V Q (T) ) 1//3 was virtually unchanged by the 
addition of P*(T). Values of the correlation coefficient Indicating the - 
quality of least-mean-square linear fit, varied from 0.99995 for the original 
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experimental data to 0.99992 with the largest added P*(T). However, the 

effect of the added P*(T) was found In changes of corresponding V , B , 

o o 

and (3B/3P) 0 with temperature. 

Thus we have a universal form for the EOS of all classes of solids 
Including the effects of temperature. In addition, there Is the capability 
(Ref. 4) to predict EOS, requiring only this single form and Input equilibrium 
data: V q , B q , and (aB/aP) Q at each temperature. While this Is a 

relatively small amount of Information necessary to predict the EOS, In the 
next section we will derive a predictive scheme which requires zero pressure 
data only at a single (reference) temperature. 

III. HIGH TEMPERATURE PREDICTIONS FROM REFERENCE TEMPERATURE 

ZERO PRESSURE DATA 

The form of the universal EOS allows us to make predictions given V , 

o 

B q , and (3B/aP) Q . Here we would like to show that for some solids, high 
temperature thermodynamical properties can be predicted via the universal EOS 
given zero pressure data at a single (reference) temperature. (Note that for 
metals, the universal binding energy relation, Eq. (2.1), provided a means of 
predicting the thermal expansion coefficient at a single temperature and zero 
pressure melting temperatures (Ref. 19)). 

The universal form (see discussion following Eq. (2.5) as well as 
Eq. (4.1)), depends at high temperatures on P*(T,V) having a much weaker 
dependence on V than - dE(V)/dV has. This characteristic seems to be 
generally true (see, e.g., Refs. 5 to 8, and 17). The new assumptions In this 
paper are that P*(T,V) Is Independent of V and linearly dependent on T 
for T >j e Q (where e Q Is the Debye temperature). 

. Experimentally, P*(T,V) can be very often assumed to be Independent of 
V. Examples are the alkali metals (Ref. 6), solid xenon (Refs. 7 and 17), 

NaCl, L1F and MgO (see Ref. 8 for the last three). However, for solid argon 


6 



and krypton the volume dependence of P*(T,V) was found (Ref. 7) to be not 
negligible, but still much weaker than that of - dE(V)/dV. Thus where 
appropriate one might assume that 


P* (T.V) -» P*(T) (3.1) 


In addition, it has been observed experimentally that above the Debye 
temperature, P*(T,V) is a linear function of T (see, e.g.. Refs. 6 to 8 
and 17). By definition. 


aP(T,V) 

3T 


3P*( T. V) > 
3T 


= * 0 (T,V) B o (T ' V) 


(3.2) 


where a o (T,V) is the volumic thermal expansion coefficient. Thus (for 

T^ e 0 ) linearity with T and Eq. (3.1) imply 

P*(T) = P*(T r ) + a d (T R ) B 0 (T r ) (T - T R ) (3.3) 

where T R > and T R Isa reference temperature. 

Combining Eqs. (2.4) and (3. 3) , we have for T ^ e D , 

P(T,V) = P(T R , V) + ol q (T r ) B o (T r ) (T - T r ) (3.4) 

By Eq. (3.4), the high temperature EOS and corresponding thermodynamical 

properties can be predicted from a knowledge of zero pressure properties at 

the reference temperature T . Specifically, P(T ,V) is specified by 

R K 

Bq(Tr), V 0 (T r ) and (3B/3P) Q (T R ) because of the universal form (Ref. 4) of EOS. 
The only additional input number is (T_). 

0 K 

From Eq. (3.4) we find that the Isothermal bulk modulus as given by 


and 


B(T,V) = B( V ) = -V 


dP(T R ,V) 

dV 


(3.5) 



V dB( V) 
" B(V) dV 


(3.6) 


The zero pressure thermal expansion can also be predicted from the EOS, 
Eq. (3.4), as 


7 



P( T r,V) ♦ a 0 (T R ) B 0 (T r ) (T - V - o (3.7) 

We will now test the predictions from Eqs. (3.4 to 3.7), making use 

(Ref. 4) of the universal EOS for the reference Isotherm, P(T R ,V). 

IV. ANALYTIC APPROXIMATION TO UNIVERSAL EQUATION OF STATE 

The universal EOS has been described In Ref. 4 and Section II. It was 

pointed out there that plots of In H(V,T) (Eq. (2.5)) versus (1 - 
1 /3 

( V/V ( T) ) ) were approximately linear at each temperature, with slope 

and Intercept yielding B q ( T) and (3B/3P) Q (T). Thus a simple, analytic 
approximation (Ref. 4) to the universal EOS Is: 


P(T,X) = 


3B.(T) 


(1 - X) exp[n (T) (1 - X)] 


(4.1) 


where 


\l/3 


X = 




(4.2) 


and 


V T > - f 


3Er 

. aP , 


(T) - 1 


(4.3) 


This analytic form, Eq. (4.1), Is particularly convenient when used as 
the reference Isotherm since It allows one to make analytic predictions of 
high temperature thermodynamical properties through Eqs. (3.4) to (3.7). In 
terms of Eq. (4.1), Eqs. (3.4) to (3.7) become, respectively, 

3B.(T ) 

P(T,X) = - s - f- (1 - X) exp[n 0 (T R ) (1 - X)] + a„(T D ) B o (T„) (T - T 0 ) 


o' R “o' 'R' 


(4.4) 


B(T,X) = 


B o<V 


[2 + (n Q (T R ) - 1) X - n 0 (T R ) X‘] exp[n 0 (T R ) (1 - X)] 


(4.5) 
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and 



IF (T.X) = 


4 + (3n 0 (T R ) - D X + n Q (T R ) (n Q (T R ) - 1 ) X 2 - n*(T R ) X 3 
3(2 + Cn 0 (T R ) - 1) X - n 0 (T R ) X 2 ) 


where, for Eqs, (4.4) to -,(4.6) . 


VV 


It Is appropriate now to compare Eqs. (4.1) and (4.4). As discussed in 
Section III, Eq. (4.4) is an approximation to Eq. (4.1) which is reasonable 
for T > 0 ^ and for those solids for which Eq. (3.1) is a good 
approximation. The numerical test described in the next to last paragraph of 


Section II suggests that Eqs. (4.1) and (4.4) are consistent with each other 


to a good approximation. As discussed earlier, the advantage of Eq. (4.4) is 
that values of B q (T), a o (T), and (3B/3 P) q (T) are required only at T = 

T R , while for Eq. (4.1) V q (T), B q (T) and (3B/3P) o (T) are required at 
all T of Interest. 


V. COMPARISON WITH EXPERIMENT 


We will now compare predictions of Eqs. (3.4) to (3.7) with experimental 
data. The analytic approximation to the universal EOS will be used for the 

reference Isotherm, Eqs. (4.4) to (4.6). Comparison will also be made, where 

.. . \ 

appropriate, with use of the Murnaghan (Ref. 20) EOS for the reference 
Isotherm in Eqs. (3.4) to (3.7). The Murnaghan EOS, which is extensively used 
(see, e.g., Refs. 21 and 22), has decreasing accuracy (Refs. 4 and 23) as 
compression Increases. 

As the universal EOS applies to all classes of solids, we carried out 
this test for representatives from three classes of solids. These are .gold, 
sodium chloride, and xenon, exhibiting respectively metallic, ionic, and 
rare-gas bonding. 
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A. Input Parameters 

Equations (4.4) to (4.6) require only three Input parameters for each 
solid. These are listed In Table I for the three representative solids, along 
with selected reference and Debye temperatures. 

A few comments must be made: for sodium chloride and xenon, the 

reference temperatures have been chosen at temperatures slightly lower than 
the Debye temperatures. This was necessary due to the very limited EOS 
experimental data, and Is not expected to Introduce Important errors In the 
predictions of the model. For gold and sodium chloride the values of 
EJ(T d ) and (3B/3P) ( T n ) have been chosen as the average of values 
taken from the references Indicated In Table I. 

It Is difficult to estimate (3B/3P) (T ) for gold from the 

0 K 

literature. Consequently, we will examine the predictions of our theory for 
two different values: the lowest one represents most of the static pressure 

experiments ( ( 3B/3P) _ ( T D ) = 5.5) and the highest value represents most of 
the ultrasonic measurements ( ( 3B/3P) (T„) = 6.5) (Ref. 24). Values for 
xenon have been taken from the analysis carried out by Birch (Ref. 25) on data 
obtained by M.S. Anderson and C.A. Swenson (Ref. 7). Finally, thermal 
expansion data have been taken from the complete compilations In Refs. 26 
and 27 for gold and sodium chloride respectively. The value for xenon has 
been taken from the study of Manzhelll et al. (Ref. 28). 

B. Prediction of High Temperature Isotherms 

In this section we will estimate high temperature Isotherms near the 
melting point from Input data measured for T « e^. The results are then 
compared with high temperature Isotherms taken from the literature. Figure 1 
compares the results for gold with the estimates of Heinz and Jeanloz 
(Ref. 24). Their estimates are not direct experimental data but 
extrapolations of the considered Isotherms using a Mle-Grunel sen equation. 
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Heinz and Jeanloz estimate that their results contain a 1 to 2 percent error 
In the pressure. A problem arises from uncertainties In the value of 
(3B/3P) (T ). In their analysis of experimental data, Heinz and Jeanloz 

0 K 

find a value close to 5.5: however, most ultrasonic measurements give values 

around 6.5 (Refs. 24, and 29 to 31). It Is Interesting to note that a fitting 

of the room temperature EOS proposed In Ref. 24 to our universal EOS (Eq. 

(4.4)) would give a value of 5.8. In the present analysis, we show the 

results for the two extreme values of 5.5 and 6.5 for (3B/3P) (T ) for 

o R 

gold (Fig. 1). The agreement Is good over the entire range of compression and 
Is considerably better than that obtained using the Murnaghan (Refs. 20 
and 22) EOS for the reference Isotherm In Eqs. (3.4) to (3.7). 

In Fig. 2 we compare the results for NaCl with the experimental data of 

Boehler and Kennedy (Ref. 32) and with data representative of the Decker model 
(Ref. 33). Again, the prediction via Eq. (4.4) Is relatively accurate. The 
prediction based on the Murnaghan reference Isotherm Is again accurate at low 
pressures but becomes Increasingly Inaccurate as the pressure Increases as Is 
to be expected. 

Finally, solid xenon, as shown In Fig. 3, confirms the results of both 
Figs. 1 and 2. Again the predictions based on Eq. (4.4) give good agreement 
when compared to experimental data (M.S. Anderson and C.A. Swenson (Ref. 7)), 
while that based on the Murnaghan P(T ,V) deviate significantly at high 

K 

compression. 

For the three different materials Investigated so far, we can conclude 
that the use of the universal description of the reference Isothermal EOS 
combined with the assumptions of Section III gives an accurate prediction of 
the high temperature EOS. We have also demonstrated that the Murnaghan EOS 
cannot describe the reference Isotherm well enough to allow extrapolations to 
high temperatures when substantial compression occurs. 
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C. Prediction of Thermal Expansion 

In Section III, a method was described for obtaining the thermal 
expansion curve given the reference temperature Isotherm (Eq. (3.7)). In this 
section, thermal expansion predictions will be compared with experiment. For 
simplicity, we will limit our study to zero pressure thermal expansion 
(Eq. (3.7)) for which experimental data are easily available In the 
literature. However, we note that the same approach predicts a at high 
pressure (P(V,T) * 0 In Eq. (3.4)), where data are generally not available. 

We will use the universal EOS for the reference isotherm, so that from 
Eq- (4.4) 

^2 (1 - X) exp[r, Q (T R ) (1 - X)] + a Q (T R ) (T - T R ) = 0 (5.1) 

X 

Here X Is defined by Eq. (4.7) and n (T ) by Eq. (4.3) evaluated at 

O K 

T r . Thus the zero pressure thermal expansion Is determined entirely by 
W’ W’ and ( 3 B/3P) 0 (T R ). 

For cubic solids, this function of volume can be reduced to a function of 
length. Figures 4 to 6 are plots of the predicted values of thermal expansion 
along with experimental data for gold, sodium chloride, and solid xenon. The 
two Input parameters can be found In Table I. Also for comparison the linear 
expansion <x q (T r ) (T - T )/3 shown - 

A common feature of Figs. 4 to 6 Is that Eq. (5.1) gives a good 
prediction of thermal expansion for temperatures ranging from the Debye 
temperature up to the melting point. The nonlinear behavior of the thermal 
expansion Is well reproduced. This Is due primarily to the universal EOS. 

Such nonlinear dependence on temperature has been difficult to predict 
(Ref. 9). 

A second common result Is that Eq. (5.1) does not accurately predict the 
thermal expansion for temperatures less than the Debye temperature. Recall 
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that Eq. (3.3) Is valid only for T > e Q and that Eq. (3.3) was assumed 

In the derivation of (5.1). Equations (3.1) and (3.3) require that a 0 (T,V) 

B Q (T, V) be Independent of T and V. Since the thermal expansion 

coefficients goes to zero at zero temperature, a Q (T,V) B Q (T,V) also 

-» 0 as T--» 0. Thus It Is to be expected that predictions from Eq. (5.1) 

would not be accurate for low temperatures. 

In Fig. 4 we Investigated the sensitivity to the Input parameter, 

(3B/3P) (T d ). Using the two extreme values for gold of 5.5 and 6.5, the 

maximum deviation between predicted thermal expansion of gold Is only 3 

percent of . Al/L(T ) . 

K 

However, Eq. (5.1) Is slightly more sensitive to the choice of the other 

parameter, a (T ). As shown In Fig. 5, a change of about 2.5 percent In 
o R 

-4 

a o (T R ), consistent with experimental accuracy ( a 0 (Tp) = 1-2x10 
-1 -4 -1 

K In Ref. 27 and 1.17x10 K In Ref. 32 for the same temperature of . 

300 K), brings maximum variations of the predicted values of about 4 percent 

of AL/L(T ). 

K 

0. Temperature Dependence of the Bulk Modulus 
In Section III we assumed that the Isothermal bulk modulus has no 
explicit dependence on the temperature and Is a function of the volume only 
(Eq. (3.5)). In Section IV, we proposed an analytical form of this function 
B(V), based on the universal EOS (Eq. (4.5)). We can now use the predictions 
of the thermal expansions obtained previously (Eq. (5.1)) to compute the 
temperature Induced volume changes and consequently, the corresponding changes 
of Isothermal bulk modulus. This computation has been done for gold, sodium 
chloride, and solid xenon and the results are plotted In Figs. 7 to 9, 
respectively. Again, the only Input parameters are the three listed In 
Table I. 
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' The results achieved with our model agree well with available 

experimental data over the Investigated temperature range. In addition to 

these general comments. Fig. 7 shows that, for gold, the choice of an Input 

value of 6.5 for (3B/3P) (T_), which corresponds to ultrasonic 

o R 

measurements, gives the best results. This was to be expected since the 

experimental data used here to test the bulk modulus prediction (from Refs. 30 

and 36) are also obtained from ultrasonic measurements. 

* 

In the case of xenon, where low temperature experimental data are 
available (Ref. 25), Fig. 9 shows a deviation of the predicted values from the 
experimental values for temperatures close to 0 K. This was also to be 
expected since the approximations made In Section III do not adequately 
represent the low temperature behavior of solids. The slope of the curves 
B q (T) must be zero at zero temperature on general thermodynamic grounds. 

This Is confirmed by the experimental data for xenon, but Is not represented 
by our model. 

E. Temperature Dependence of (3B/3 P) q 
T he same approach used for the Isothermal bulk modulus can be used for 

(3B/3P) since this quantity Is also assumed to have no explicit 
o 

temperature dependence (Eq. (3.6)). We use Eq. (4.6) and the previously 
calculated values of thermal expansion (Eq. (5.1)) for sodium chloride and 
xenon to obtain a prediction of the thermal behavior of (3B/3P) o< 

Figures 10 and 11 present respectively the results for sodium chloride 
and xenon, obtained by using the three Input parameters listed In Table I. 
These results are compared with experimental data from ultrasonic measurements 
(sodium chloride) or from analysis of high pressure EOS data (xenon). 
Unfortunately, to our knowledge, no similar high temperature set of data 
exists for gold. 
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For sodium chloride, the agreement between the prediction of Eq. (4.6) 

and the experimental data, as presented In Fig. 10, Is good. For xenon 

(f^g* 11), agreement Is worse. However, the accuracy of the experimental data 

which Is extrapolated from static high pressure measurements Is questionable. 

For example, the two sets of data presented In Fig. 11 are extrapolated from 

the same set of high pressure experiments, using two different methods 

(Refs. 7 and 25). As a consequence, It Is clear that one must add, the 

uncertainties of the high pressure measurements to the uncertainty caused by , 

the extrapolation method used to obtain the values of (3B/3P) . Therefore, 

. o 

a check of our predictions of the thermal variations of (aB/aP) Q will 
require more complete and accurate sets of experimental data. 

VI. SUMMARY 

We have discussed a recently proposed "universal" EOS which accurately 
describes Isothermal P-V data for all classes of solids In compression. Phase 
transitions were excluded from the analysis. Results from Ref. 4 show that 
the EOS can be predicted from the universal relation If V q (T), B q (T) and 
(3B/aP) o (T) are known for the temperatures of Interest. 

In this paper we have used the fact that an additional assumption Is 
often valid. Namely It was assumed that the thermal pressure Is Independent 
of V and varies linearly with T for T > e^. With solids for which 
this Is an accurate approximation we showed that high temperature 
thermodynamical properties can be predicted from a Q (T R ), V q (T r ), 

B IT d ) and (3B/3P) (T D ) , where T„ Is a single reference 
temperature. An analytic approximation was Introduced for the universal EOS. 
This led to a number of simple formulas characterizing the high temperature 
EOS In the absence of phase transitions. Simple expressions of particular 
Interest were given for: (1) the high temperature Isotherms; (2) the 

temperature dependence of the thermal expansion for a given P; and (3) the 
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temperature dependence of the P = 0 bulk modulus and Its pressure 
derivative. These simple expressions were compared with experiment for gold, 
sodium chloride, and xenon. Good agreement between prediction and experiment 
was obtained for the high temperature Isotherms, as well as the temperature 
dependence of the thermal expansion, bulk modulus, and pressure derivative of 
the bulk modulus for gold, sodium chloride, and xenon. 

Finally, we want to emphasize the wide generality of a model which 
describes the high temperature thermodynamic behavior of a metal and an alkali 
halide as well as a rare gas solid, from P = 0 Input data at a single 
(reference) temperature. 
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TABLE I. - SUMMARY OF THE THREE INPUT PARAM- 
ETERS AS WELL AS DEBYE TEMPERATURES 


0 D AND REFERENCE 
TEMPERATURES T 

K 


Input data 

Gold 

NaCl 

Xenon 

B o (T R )(10l0 P a) 

«o(Tr)(10-5k-1) 
( 3B/3P) 0 { T r ) 

a 1 6 . 6 
d 4.25 
35.5-6.5 

t>2 . 35 

ei2.0 

b5.35 

CO. 302 
f 60 .0 

C7.8 

Temperatures 




300 

298 

6 

9 

162 

320 

6 



a Reference 24. 
b Reference 4. 
c Reference 25. 

d Reference 26, data taken at 293 K. 
e Reference 27, data taken at 293 K. 
f Reference 28. 

9 Reference 3. 



















PRESSURE, 10 8 Pa PRESSURE, 10 8 Pa 



RELATIVE VOLUME V 

V o (300 K) 


Figure 1. - Prediction of a high temperature PIT, V) 
isotherm: Gold at 1000 K. Comparison of predic- 
tions based on universal P(T R ,V)(eq. (4. 41) and on 
the Murnaghan model for P(T R , V) with experimen- 
tal data from Heinz and Jeanloz (ref. Z4I. 



V 0 (298 K) 

Figure 2. - Prediction of a high temperature P(T,V) iso- 
therm: NaCI at 773 K. Comparison of predictions based 
on universal P(T R ,V)(eq. (4. 4)) and on the Murnaghan 
model for P(T R , V) with experimental data from Boehler 
and Kennedy (ref. 32) and representative points from 
the Decker EOS (ref, 33). 



THERMAL EXPANSION AL/L(T P ),% 



RELATIVE VOLUME V 

V 0 160 K) 

Figure 3. - Prediction of a high temperature PIV) iso- 
therm: Xenon at 159 K. Comparison of oredictions 
based on universal PlVVIleq. (4.4)1 and on the 
Murnaghan model for P(T R . V) with experimental 
data from M.S. Anderson and Swenson (ref. 7). 



Figure 4. - Prediction of thermal expansion of gold: Plot 
of prediction from equation 15. 1) with smoothed experi- 
mental data from reference 26. To examine how well 
nonlinear expansion is predicted, the linear result is 
also plotted as (T - T r )(^j(T r V3. 


THERMAL EXPANSION AL/L (T R ), % THERMAL EXPANSION AL/L(T R ),% 



Figure 5. - Prediction of thermal expansion of sodium 
chloride plotted as in figure 4 with smoothed experi- 
mental data from reference 27. 



TEMPERATURE, K 


Figure 6. - Prediction of thermal expansion of solid 
xenon plotted as in figure 4 with experimental 
data from reference 34 (T < 100 K) and reference 
35 (T> 100 K). (See also ref. 7.) 


ISOTHERMAL BULK MODULUS, 10 10 Pa ISOTHERMAL BULK MODULUS, 10 JU Pa 



Figure 7. - Prediction of temperature dependence of the 
isothermal bulk modulus of gold. Plot of the present 
model legs. (4. 5) and (5. D) with experimental data 
from references 30 and 36. 
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Figure 8. - Prediction of temperature dependence of 
the isothermal bulk modulus of sodium chloride.- 
Plot of the present model (eqs. 16. 5) and (5. U) 
with experimental data from references 32 and 37. 



XENON 



Figure 9. - Prediction of temperature depen- 
dence of the isothermal bulk modulus of 
solid xenon: Plot of the present model 
legs. (4. 5) and (5. D) with experimental 
data from reference 25. 


NaCI 



Figure 10. - Prediction of temperature dependence 
of (dB 0 (T, P)/dP) T for sodium chloride: Plot of 
the present model (eqs. 14. 6) and (5. 11) with 
experimental data from references 37 and 38. 



TEMPERATURE, K 


Figure 11. - Prediction of temperature dependence of 
(dB„(T, PF/dPK for solid xenon: Plot of the present 
model leas. (4. 5) and (5. lilwith experimental data 
extrapolated from high pressure measurements in 
references 7 and 25. 
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